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ABSTRACT 

A spectroscopic analysis of the DBA (or DAB) white dwarf KUV 02196+2816 
is presented. The observed hydrogen and hehum hne profiles are shown to be 
incompatible with model spectra calculated under the assumption of a homoge- 
neous hydrogen and helium chemical composition. In contrast, an excellent fit 
to the optical spectrum of KUV 02196+2816 can be achieved if the object is 
interpreted as an unresolved double degenerate composed of a hydrogen-line DA 
star and a helium-line DB star. The atmospheric parameters obtained from the 
best fit are T^s = 27, 170 K and logg = 8.09 for the DA star, T^s = 36,340 K 
and logg = 8.09 for the DB star, which implies that the total mass of the system 
(M ~ 1.4 Mq) is very close to the Chandrasekhar limit. Moreover, the effective 
temperature of the DB stars lies well within the so-called DB gap where very 
few bright DB stars are found. The implications of this discovery with respect 
to the DAB and DBA spectral classes and to the evolution of double degenerate 
binaries are discussed. 

Subject headings: binaries: spectroscopic, stars: individual (KUV 02196+2816), 
white dwarfs 



1. Introduction 

The discovery of double degenerate binary systems is crucial to our attempt at better 
constraining their origin and evolution. In particular, these binary systems are thought to 
be the progenitors of Type la supernovae if the total mass of the two components exceeds 
the Chandrasekhar limit and if their orbital period is small enough to allow them to merge 
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within a Hubble time. However, onl y a small percen t age o f do uble degenerate binaries lead 
to supernova explosions according to lNelemans et al.l (120011 ) and lNapiwotzki et al.l (I20011). In 



2006 , only ten such binaries composed of two white dwarfs were known (Ivan der Sluys et al. 
20061 ). Since then, the SPY (ESO Supernovae la Progenitor surve Y) radial velocity surve y 



has observed more than 1000 white dwarfs and pre- white dwarfs (INapiwotzki et al.l 120071 ). 
One of the conclusions of that survey is that double degenerate binaries cannot explain the 
existence of every low mass white dwarf (0.45 Mq or less). Indeed, only 42% of low-mass 
He-core white dwarfs are found in close binary systems. 

Double degenerate binaries are difficult to detect from single slit spectroscopy alone. 
This is particularly true for system composed of two DA stars since the composite spec- 
trum can be reprod uced almost perfectly with a single hydrogen-rich model atmosphere 
(ILiebert et al.lll99ll). In order to detect such spectroscopically invisible DA+DA systems, 
Lajoie fc Bergeron! (120071 ) compared effective temperatures determined from optical and ul- 
traviolet (lUE) spectra, and identified three candidate systems. Double degenerate systems 
composed of white dwarfs with different spectral types are easier to detect, however, since 
the modeling of these objects under the as s umpt ion of a single star leads to spurious results. 
This is precisely how [Bergeron fc LiebertI (120021 ) showed that the DAB PG 1115+166 was 
in fact a double degenerate binary system composed of a DA and a DB star. Their analysis 
revealed that the optical spectrum could not be reproduced with a single model atmosphere 
with a mixed hydrogen and helium composition, or even a stratified chemical composition. 
Instead, the authors showed that a combination of a DA and a DB model spectrum could 
perfectly match the observed line profiles and photometry. Additional objects are also dis- 



the PG catalog (iGreen et al. 



cussed in iPereira et al.l (120051) and references therei n. PG 1115-H66 wa s first identified in 



19861 ) as a DA5. Later, iMaxted et al.l (120021 1 confirmed the bi 



nary nature of PG 1115+166 and obtained an orbital period of 30.09 da ys, too long to allow 



the w hite dwarfs to merge within a Hubble time. The binary models of Ivan der Sluys et al. 



(120061 ) explained that the older nature of the DB star in PG 1115+166 is the result of the 
loss of its hydrogen-rich outer layers after it has evolved into a giant phase. 

We are currently conducting a spectroscopic survey of all the white dwarfs discovered 
in the Kiso Ultraviolet (KUV) sur vey with the airn of red etermining the luminosity function 
using the spectroscopic method of iBergeron et al.l (Il992l . see also Liebert et al. 2005). One 
of these white dwarfs, KUV 02196+2816 (WD 2 19+2 82, V = 17.30) was spectroscopically 
identified as a DBA star by iDarling &: Wegnen (119961 ). However, this object should have 
been classified as a DAB star since the hydrogen lines are actually stronger than the helium 
lines (see § 2 below). DAB white dwarfs are particularly interesting since they may represent 
the key to our understanding of the evolution of DAO and hot DB sta rs, and of the existe nce 
of a strong deficiency of DB stars between 30,000 and 45,000 K (see lVennes et al.ll2004l . for 
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a review 



We report in this paper the discovery that KUV 02196+2816 is actually a double degen- 
er ate binary composed of a D A and a DB star. An analysis identical to that of PG 1115+166 



by [Bergeron fc LiebertI (120021 ) indeed reveals that the spectrum of KUV 02196+2816 cannot 
be reproduced by a single star model, and that a composite DA and DB model spectrum 
provides a much better fit to the observations. We first present our spectroscopic observa- 
tions in § 2. The model atmosphere calculations described in § 3 are then used in § 4 to 
analyze in detail the optical spectrum of KUV 02196+2816 using homogeneous as well as 
composite model atmospheres. Our discussion follows in § 5. 



2. Spectroscopic Observations 



Optical spectroscopy for KUV 02196+2816 has been obtained on 2007 December 6 using 
the Steward Observatory 2.3-m telescope equipped with the Boiler & Chivens spectrograph 
and a Loral CCD detector. The 4. "5 slit together with the 600 line mm~^ grating in first 
order provided a spectral coverage of 3200-5300 A at an intermediate resolution of ~ 6 A 
FWHM. 

Our optical spectrum for KUV 02196+2816 is displayed in Figure [TJ The hydrogen 
Balmer lines are clearly stronger than the helium lines and the star should thus be classified 
as a DAB white dwarf. The s pectrum of KUV 02 196+2816 is also compared with CD 323, th e 



as a JJAn wnite awari. i ne s pectrum oi KU v u/ i^D+zaiD is also comparea witn KjU o/j, tn ( 
prototype of the DAB c lass jpereira et allboosh . PG 1115+166 feergeron fc LiebertI [2002[ ) 
and MCT 0453-2933 JWesemael et al.lll994h . PG 1115+166 has already been discussed 
above while MCT 0453—2933 has been analyzed by IWesemael et al.l (119941 ) and was shown 
to be an unresolved composite system consisting of a DA white dwarf together with a DB 
or a DBA star. F r om re peated spectroscopic observations at Ha as part of the SPY survey. 



Napiwotzki et al.l (120051 ) found variations of the line core shap es. In one of the spectra, they 
noticed a splitting at Ha (see also Fig. 7 of IVoss et al.l 120071 ). indicating that both stars in 



the binary system contain considerable amounts of hydrogen. It was concluded that MCT 
0453—2933 is composed of a DA and a DBA, or even a DAB star. 



m 



The similarity between M CT 0453—2933 and PG 1115+166 has already been discussed 



Bergeron fc LiebertI (120021 ). where the double degenerate binary nature of PG 1115+166 



was unveiled. The spectrum of KUV 02196+2816 shows the same He I features as these 
two stars: the 4026 and 4471 A absorption lines are about the same strength as those of 
PG 1115+166 while the features at 4713, 4921 and 5015 A are only slightly stronger. The 
hydrogen Balmer lines are definitely weaker in the spectrum of KUV 02196+2816 than in 
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Fig. 1. — Comparison of our optical spectrum of the DAB star KUV 02196+2816 with those 
of the DAB stars G D 323, PG 1115+166 and MC T 0453-2933. The last two objects have 
been interpreted by lBergeron fc LiebertI (120021 ) and lWesemael et al.l (119941 ). respectively, as 
unresolved composite systems consisting of a DA white dwarf and a DB or DBA companion. 
The spectra are normalized at 4400 A and are shifted vertically by 0.5 for clarity. 



PG 1115+166, however, and more typical of the strength observed in GD 323. 



3. Model Atmospheres and Synthetic Spectra 



The model atm ospheres and synthetic spectra for DA stars are described at length in 



Liebert et al.l (120051) and references therein. The models for DB and DBA stars used in this 



analysis rely on a modified version of our model atmosphere code for DA stars in which we 

lelium opacity sources from the DB model atmosphere code described 



have included all the 
in 



Beauchamp et al 



of iBeauchamp et al. 



1996r) . and in particular the improved Stark profiles of neutral helium 



(119971 ). The mixed hydrogen and helium models assume a homogeneous 
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chemi cal composition. We refrain here from using stratified models since iBergeron &: Liebert 
( 120021 ) demonstrated that such models could not improve the fit to the optical spectrum 
of PG 1115+166 with respect to homogeneous models. Also, because the hydrogen and 
helium lines in the spectrum of KUV 02196+2816 are almost of equal strength, DBA models 
with unusually high helium abundances had to be calculated for this analysis. Our DBA 
model grid covers a range of Teg = 14,000 (2000) 50,000 K, log^ = 7.0 (0.5) 8.5, and 
log A^(He)/iV(H) = —1.0 (0.5) 1.0, where the numbers in parentheses indicate the step value. 
Pure helium models were calculated as well. 



4. Model Atmosphere Analysis 



Our fitting te chnique relies on the nonlinear least-squares method of Levenberg-Marquardt 
(jPress et al.lll986l ). which is based on a steepest descent method. The model spectra (con- 
volved with a Gaussian instrumental profile) and the optical spectrum of KUV 02196+2816 
are first normalized to a continuum set to unity. The calculation of is then carried out 
in terms of these normalized line profiles only. All atmospheric parameters - T^g, logg and 
iV(He)/A^(H) - are considered free parameters. When fitting DA+DB model spectra, the to- 
tal flux of the system is obtained from the sum of the monochromatic Eddington fluxes of the 
individual components, weighted by their respective r adius. The stel l ar rad ii are obtained 
from evolutionary models similar to those described in lFontaine et al.l (120011 ) but with C/0 
cores, g(He) = logMne/M* = 10~^ and g(H) = 10""^, which are representative of hydrogen- 
atmosphere white dwarfs, and qCRe) = 10~^ and g'(H) = 10~^°, which are representative of 
helium-atmosphere white dwarf^j. 

Our best flts to the optical spectrum of KUV 02196+2816 using homogeneous and 
composite DA+DB models are shown in Figure El The solution using h omogeneous models 
depicts the same problem as that mentioned in iBurleigh et al.l (120011 ) with the flt to PG 
1115+166: the He ll A4686 feature predicted by our models is not observ ed in the spectrum 
of KUV 02196+2816. We note, however, that iBergeron &: LiebertI (120021 ) did not encounter 
this problem while fltting the spectrum of PG 1115+166, even though this particular line 
was included in their model calculations. Our homogeneous solution for KUV 02196+2816 is 
~ 3000 K hotter and ~ 0.2 dex lower in logg than that obtained for PG 1115+166, but more 
importantly, the helium abundance is about a factor of 40 larger, and our solution predicts 
a weak but detectable He ll absorption feature that is simply not observed. We also see in 
Figure [2] that the low hydrogen Balmer lines (H/9 to B.6) are predicted too strong while the 



http://www.astro.uinontreal.ca/~bergeron/CoolingModels/ 
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higher members are too weak. The homogeneous model yields a satisfactory fit for the He I 
lines AA4026, 4713, and 5015, but a poor fit to the remaining helium lines. 



KUV 02196+2816 




4000 4200 4400 4600 4800 5000 

Wavelength (A) 

Fig. 2. — Our best fits to the optical spectrum of KUV 02196+2816 using homogeneous and 
composite DA+DB models. The atmospheric parameters for each solution are given in the 
figure. Both the observed and theoretical spectra are normalized to a continuum set to unity. 
The spectra are shifted by a factor of 0.4 from each other for clarity. Clearly, the DA+DB 
solution provides the best fit to the overall spectrum. 



Our grid of model atmospheres with mi xed H/He compositions has already been ap- 
plied successfully to fit DB and DBA stars (iBeauchamp et al.lll996l ) or hotter DAO stars 
( jBergeron et al.lll994l ). Hence the bad fit displayed at the t op of Figure [2] do e s not refiect the 
inability of our models to fit this star. We also note that lEisenstein et al.l (120061 ) have ob- 
tained good fits to DB stars found in the Sloan Digital Sky Survey (SDSS) with temperatures 
similar to that of KUV 02196+2816 {T^s = 35, 200 K), although none of these hot DB stars 
had hydrogen abundances as high as that inferred here (A^(H)/A^(He) ~ 0.7). Even if we ig- 
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Fig. 3. — Relative energy distributions for our best composite DA+DB fit displayed in 
Figure El The thin lines show the individual contributions of the DA and DB components, 
properly weighted by their respective radius, while the thick line corresponds to the total 
monochromatic flux of the composite system. 



nore the bad quality of the fit depicted at the top of Figure El the atmospheric parameters we 
obtain for KUV 02196+2816 based on homogeneous models represent a problem from an as 



trophy sical point of view. Indeed, the typical hydrogen abundances measured by IVoss et al. 



(120071 ) in DBA stars from the SPY survey are of the order of A^(H) /A^(He) ~ 10"^ to 10"^ (see 
their Table 2) — or even 10~^ in some extreme cases — , and most of them have effective 
temperatures well below 30,000 K. Similarly, the helium abundances measu red in typical 



DAP and DAB stars range from A^(He)/A^(H) ~ 10"^ to 0.1 (see Fig. 7 of IVennes et al. 



20041 ). It is thus difficult to reconcile the atmospheric composition determined here for KUV 



02196+2816 using homogeneous models with other single DAB/DAO or DBA white dwarfs. 

In contrast with the homogeneous solution, our DA+DB solution shown at the bottom 
of Figure El provides an excellent fit to the Balmer and neutral helium lines simultaneously. 
Our DA+DB fit yields a value of 0.543, a value that is significantly lower than that 
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obtained for the homogeneous solution, = 1.416. All observed features are reproduced 
in detail and this is clearly a much better solution for KUV 02196+2816. The effective 
temperature determined for the DA component is significantly lower than the value achieved 
under the assumption of a single object with a homogeneous composition. This is a direct 
consequence of the hydrogen lines being diluted by the continuum flux of the DB star (see 
below); the Balmer lines are weakened, and the effective temperature of the model needs 
to be increased to match the observed line profiles. The effective temperature of the DB 
component. Teg = 36, 340 K, is higher than the value obtained with our homogeneous models, 
Teff = 35, 200 K, yet the He ll A4686 feature only appears as a small depression in the blue 
wing of He I A4713. Again in this case, the presence of the DA star dilutes all the helium 
lines present in the DB spectrum. But the main physical reason for this difference in strength 
is that the overall opacity of the mixed H/He model is lower than that of the pure helium 
model at the same effective temperature, resulting in lower atmospheric pressures, which in 
turn favors the ionization of helium in the mixed atmosphere. Hence the He ll A4686 feature 
appears stronger in the mixed model than in the pure helium model, even though the former 
is ~ 1000 K cooler. 

The surface gravities of both components of the system are almost identical, logg = 
8.089 for the DA and logg = 8.095 for the DB. These values can be translated into masses 
and cooling ages using the evolutionary models described above. We obtain respectively for 
the DA and the DB star 0.69 and 0.68 M©, and cooling ages of 1.76 x 10"^ and 5.13 x 10^ 
years. Since both stars have the same radius, the contribution of each component to the total 
luminosity is only a function of the effective temperature. Since the DB star is significantly 
hotter than the DA component, the former will contribute more to the combined luminosity 
of the system. This is illustrated quantitatively in Figure [3] where the contribution of each 
component to the total flux is depicted. There is a significant contribution of the flux of the 
DB star in the optical regions of the energy distribution. In particular, the cores of the lower 
Balmer lines are filled in by the continuum flux of the DB star, resulting in the poor fits of 
the homogeneous solution displayed in Figure [2] and discussed above. 



5. Discussion 

Our analysis has shown that the simultaneous presence of hydrogen and helium in the 
spectrum of KUV 02196+2816 is better explained in terms of an unresolved binary system 
composed of a DA white dwarf and a DB star. We note that in at least two instances, 
PG 1115+166 and MOT 0453—2933 (displayed in Fig. [1] and discussed above), this binary 
interpretation has been confirmed through spectral velocity variations, hence similar obser- 
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vations of KUV 02196+2816 should also confirm its binary nature. Despite the evidence of 
our spectroscopic analysis, we discuss in the following the difficulty with interpreting KUV 
02196+2816 as a single star, even from an evolutionary point of view. 

There are actually two scenarios that can produce mixed hydrogen and helium at- 
mospheres in the temperature range considered here. The first scenario, discussed in the 
context of DBA stars, is the accretion of hydrogen from the interstellar medium onto a 



helium -dominated atmosphere. However, the stellar wind model proposed by lFontaine et al. 
( I2OO5I ) might prevent the accretion of hydrogen for stars hotter than T^s ~ 20, 000 K. Hence 



the presence of hydrogen in KUV 02196+2816 cannot be explained by this scenario, and 
another explanation must be sought for the existence of KUV 02196+2816 if it is to be 
interpreted as a single star. 



Another scenario has been proposed by iFontaine fc Wesemaell (119871 ). where residual 
amounts of hydrogen left in the envelope of hot DO stars would gradually diffuse to the 
surface as the white dwarf evolves along the cooling sequence. This build up of hydrogen at 
the photosphere would gradually turn DO white dwarfs into DA stars. Prior to the SDSS, 
the absence of any helium- rich atmosphere white dwarfs below T^s ~ 45, 000 K (and above 
30,000 K) would impose a limit of Mh ^ 10^^^ Mq to the total mass of hydrogen left after 
the post-AGB phase. Such minute amounts of hydrogen are actually sufficient to turn a// hot, 
helium-rich atmosphere white dwarfs into DA stars by the time they reach T^s ~ 45, 000 K, 
which defines the blue edge of the so-called "DB gap", a range in effective temperature 
between T^g ~ 30, 000 and 45,0 00 K where 110 heli um atmosphere white dwarf (DO or DB 
stars) had ever been identified (ILiebert et al.lll986l ). 



More r ecently, however, the SDSS has revealed the existence of several hot DB stars 
in this gap (lEisenstein et al.ll2006l ). although the fraction of helium-dominated atmospheres 
in this temperature range remains significantly lower than that found at higher or lower 
temperatures . We n ote that SDSS white dwarfs are much fainter than those investigated by 
Liebert et al.l (119861 ) in the Palomar-Green sample. The implication of this result is that a 
few white dwarfs must necessarily survive the DO to DA transition at high temperatures. 
This in turn implies that the amount of hydrogen left in the envelope of pre-white dwarfs 
during the post-AGB phase is even smaller than previ ously anticipated. Wemust conclude 
that the born-again post-AGB scenario proposed by IWerner fc Herwigl (120061 ) in which a 
violent mixing event is induced by a late helium fiash in the post-A GB phase is e f ficient 
enough to leave virtually no hydrogen behind. The recent discovery by iDufour et al.l (120071 . 
20081 ) of a new class of white dwarf stars with carbon-rich atmospheres between ~ 18, 000 
and 24,000 K even suggests that these fiash events might be efficient enough to leave no 
helium behind! In this context, it is perhaps no longer surprising to find hot DB stars in 
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the gap. These hot DB stars are most hkely the immediate progenitors of the carbon-rich 
atmosphere white dwarfs discovered by Dufour et aL 

The significant increase in the number of hehum-atmosphere DB white dwarfs below 
Tgff ~ 30, 000 K — whether in the PG or SDSS samples — can only be explained in terms 
of the convective dilution of the superficial hydrogen atmosphere by the underlying helium 
convective envelope, provide d that the hydrog en layer is sufficiently thin {M\\lMtnt ~ 10^^^, 
MacDonald k Venneslll99ir i. As discussed bv Bergeron fc LiebertI (|2002|), if the DO-to-DA 
and DA-to-DB transition scenarios discussed above are correct, white dwarfs with mixed 
hydrogen and helium compositions must be sought either at the hot end of the DB gap (or 
deficit) near 45,000 K where an extremely thin hydrogen atmosphere in diffusive equilibrium 
on top of the helium envelope would make the star appear as a DAO star, or near the cool 
end of the gap, where convective dilution of the thin hydrogen atmosphere occurs. Excellent 
candidates of such stars caught in the act of changing from one spectral type to another are 
PG 1305—017, a DAO star at Tpff = 44, 00 K whose spectrum is better reproduced with 
stratified atmosp heres (IBergeron et al.lll994j ). and GD 323 (shown in Fig. 1), a DAB star at 
Teff = 28, 750 K (jKoester et al.lll994j ) whose spectrum exhibits spectroscopic variations that 
have been interpreted as surface abundance inhomo geneities resulting from the convective 
dilution of hydrogen into helium (jPereira et al.ll2005l ). 



Going back to KUV 02196+2816, its effective temperature near 35,000 K obtained under 
the assumption of homogeneous single-star models does not fit well into this picture. First, 
the convective efficiency of the helium envelope at that temperature (the He ll convection 
zone in this case) is much too low to produce any dilution of the hydrogen superficial atmo- 
sphere. And second, since KUV 02196-1-2816 is a full 10,000 K cooler than the blue edge of 
the DB gap near 45,000 K, it has necessarily survived the DO-to-DA transition. This in turn 
implies that the total amount of hydrogen left in the envelope during the post-AGB phase 
must be extremely small, and certainly too small to account for the hydrogen abundance 
inferred from our homogeneous model atmosphere analysis. 

Thus, the DA+DB solution proposed here for KUV 02196-f2816 not only provides a 
much better fit to the optical spectrum than the single star model, but it also represents the 
only viable solution from an evolutionary point of view. Interestingly enough, the effective 
temperature we infer for the DB component of the system, T^s = 36, 340 K, puts it right in 
the middle of the DB gap. This makes the DB component of the KUV 02196+2816 system 
the brightest DB white dwarf ever identified in the gap. From the atmospheric parameters of 
the DA and DB stars, we estimate the absolute visual magnitude of the system at My = 9.51, 
which combined with the visual magnitude of V = 17.30 yields a distance of 360 pc. 



Wachter et al 



(120031 ) also found from the 2MASS photometry that KUV 02196+2816 



was a good binary candidate compose d of a white dwarf a nd a low-mass main-sequence star. 
This binary nature was confirmed by iFarihi et al.l (120061 ) in their study of white dwarf-red 
dwarf systems resolved with HST. Farihi et al. indeed found that KUV 02196-1-2816 has a 
red dwarf companion and both stars are believed to be gravitationally bound. Since KUV 
02196-1-2816 itself is a double degenerate binary, the system is thus composed of three stars, 
two of which are unresolved . As mentioned above, h igh-resolution spectroscopic observations 
similar to that obtained by iNapiwotzki et al.l (120051 ) for the double degenerate system MCT 
0453—2933 could help confirm our binary interpretation. Finally, we note that since the 
total mass of the KUV 02196+2816 system (M - 1.4 M©) is close to the Chandrasekhar 
limit, a precise measurement of its orbital period would establish whether KUV 02196-1-2816 
represents a likely supernova candidate. 
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